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Abstract 

We estimate the electromagnetic decay rates of XcO ~~ * 77 an d Xc2 ~* 77 by taking into 
account both relativistic and QCD radiative corrections. The decay rates are derived in the 
Bethe-Salpeter formalism and the QCD radiative corrections are included in accordance with 
factorization assumption. Using a QCD-inspired interquark potential, we obtain relativistic BS 
wavefunctions of XcO and \c2 by solving BS equations for the corresponding 2S+1 Lj states. Our 
numerical result for the ratio R = ^[ XcQ ~ >77 | is about 11 — 13 which agrees with the update 

E760 experimental data. Explicit calculations show that the relativistic corrections due to spin- 
dependent interquark forces induced by gluon exchange enhance the ratio R substantially and 
its value is insensitive to the choice of parameters that characterize the interquark potential. 
Our expressions for the decay widths are identical with that obtained in the NRQCD theory to 
the next-to-leading order in v 2 and a s . Moreover, we have determined two new coefficents in 
the nonpcrturbative matrix elements for these decay widths. 
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1 Introduction 



Charmonium physics is in the boundary domain between perturbative and non-perturbative 
QCD. Charmonium decays may provide useful information on understanding the nature of in- 
terquark forces and decay mechanisms. Both QCD radiative corrections and relativistic correc- 
tions are important for charmonium decays, because for charmonium the strong coupling constant 



a s ~ 0.3 [defined in the MS scheme (the modified minimal subtraction scheme)] and the velocity 
squared of the quark in the meson rest frame v 2 ~ 0.3, both are not small. Decay rates of heavy 
quarkonium in the nonrelativistic limit with QCD radiative corrections have been studied (see, 
e.g.,ref.Jl|, [2], ||[ |[]). However, the decay rates of many processes are subject to substantial relativis- 
tic corrections Q. With this goal in mind, people have studied relativistic corrections to the decay 
rates of S-wave charmonium 7] c , J/ip and their radial excited states |5|, [6|, 0]. These results show 
that relativistic effects are significant in the cc systems especially for the hadronic decays of J/ip. 
In the present paper, we will investigate the relativistic corrections to the electromagnetic decays 
of P-wave charmonium states XcO 77 and Xc2 — ► 77- 

The P-wave charmonium decays are interesting. Now their experimental results are quite 



uncertain. The Crystal Ball group (see j|, 16] and references therein) gives T(xc0 — ► 77) = 4.0 ± 
2.8keV. But for T(x C 2 - > 77) , its central value differs significantly among various experiments 



1C, 11 1, and the ratio of the photonic widths for XcO an dXc2 states measured by E760 is much 



larger than that measured by other two groups. Theoretically, in the nonrelativistic limit, the ratio 
is T 01' anc ^ *t will increase to about 7.4 [|| if QCD radiative corrections are considered. Recently 
a rigorous factorization formula which is based on NRQCD has been developed for calculations of 
inclusive decay rates of heavy quarkonium. In this approach the decay widths factor into a set of long 
distance matrix elements of NRQCD with each multiplied by a short distance coefficient. To any 
given order of relative velocity v of heavy quark and antiquark, the decay rates are determined by 
several nonperturbative factors which can be evaluated using QCD lattice calculations or extracted 
by fitting the data. The study of the photonic decays of XcO and Xc2 can also povide a determination 
for the nonperturbative factors in the decays of P-wave quarkonium. 



In this paper, we will use the Bethe-Salpeter (BS) formalism [12| to derive the decay am- 
plitudes and to calculate the decay widths of XcO — ► 77 and Xc2 —* 77- The meson will be treated 
as a bound state consist of a pair of constitutent quark and antiquark (higher Fock states such as 
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\QQg > and \QQgg > are neglected because they don't contribute to electromagnetic decays) and 
described by BS wavefunction which satisfies the BS equation. A phenomenological QCD-inspired 
interquark potential will be used to solve for the wavefunction and to calculate the decay widths. 
Both relativistic and QCD radiative corrections to next-to-leading order will be considered based on 
the factorization assumption for the long distance and short distance effects. The remainer of this 
paper is organized as follows. In Sec. 2 we derive the reduced BS equation for any angular momen- 
tum state 2S+1 Lj of heavy mesons. In Sec. 3 we give out the decay amplitudes of \J ~^ ll{J = 0) 2) 
and use the solved relativistic BS wavefunctions to calculate the numerical results of decay widths. 
A summary and discussion will be given in the last section. 

2 Reduced BS equations for any angular-momentum state 2S+1 Lj 
of heavy mesons 

Define the Bethe-Salpeter wavefunction, in general, for a QiQ 2 bound state IP > with overall 



mass M and momentum P = (yP 2 + M 2 , P) 

x ( Xl ,x 2 ) =< 0\TMx 1 )Mx 2 )\P >, 
and transform it into momentum space 

XP (p) = e- iPX ( d 4 xe-^ x X (xi,x 2 ). 



Here pi{m\) and p 2 (m 2 ) represent the momenta(masses) of quark and antiquark respectively, 

X = 7/1X1 + 772^2, x = xi-x 2 , 

P =Pl+P2, P= V2P1 ~ VlP2, 

where 7^ = ^^(7 = 1,2). 



We begin with the bound state BS equation [12| in momentum space 



- mi)xp(p)&fc + m 2 ) = — J <fkG(P,p- k) X p(k), (3) 

where G(P,p — k) is the interaction kernel which dominates the interquark dynamics. In solving 
(H), we will employ the instantaneous approximation since for heavy quarks the interaction is 



dominated by instantaneous potentials. Meanwile, we will neglect negative energy projectors in 
the quark propagators which are of even higher orders. Defining three dimensional BS wavefunction 

®p(P) = J dp°xp(p), 
we then get the reduced Salpeter equation for $p(p*) 

(M — Ei — E 2 )®(p) = A^ 7o J d 3 kG(P,p- £)$(£) 7o A 2 • (4) 

Here G(P,p — k) represents the instantaneous potential, A + (A_) are the positive (negative) energy 
projector operators for quark and antiquark respectively 

A i = Ei + 707 • pi + m7o 



; 2E 



2 = E 2 - 7o7 • P2 ~ mio 
2E 2 

Ei = \Jtfi +m\, E 2 = \jp 2 +m%. 

We will follow a phenomenological approach by using QCD inspired inter-quark potentials, which 
are supported by both lattice QCD and heavy quark phenomenology, as the kernel in the BS 
equation. The potentials include a long-ranged confinement potential (Lorentz scalar) and a short- 
ranged one-gluon exchange potential (Lorentz vector) 

V(r) = Vs(r)+-y lt ®'fV v (r), 
(1 - e~ ar ) 



Vs(r) = Xr- 



ar 



Vy{r) = -^e— , (5) 

where the introduction of the factor e~ ar is to regulate the infrared divergence and also to in- 
corporate the color screening effects of dynamical light quark pairs on the QQ linear confinement 
potential. In momentum space the potentials become 

G(p) = Gsift + i^i^Gvip), 
G S (P) = --S 3 (P)+ X ' 



a it 2 (p 2 + a 2 ) 2 ' 

GvW = (6) 



where ct s (p) is the quark-gluon running coupling constant and is assumed to become a constant of 
O(l) asp 2 ^ 

M = 127T 1 
QaW 27 /n(a+^/A^ OT )- 

The constants A, a, a and Aqcd are the parameters that characterize the potential. 

For any given angular-momentum state 2S+1 Lj of mesons, its three dimensional wavefunc- 
tion in the rest frame of mesons takes the following two forms: 

(i) S=0, then J=L, 

<$> Lm {p) = A+7o(l + 7o)757oA?.YL m (p)^(p); (7) 

(ii) S=l, then J=L-1,L,L+1 for L/0 or J=l for L=0, 

$Jm(p) = E < JM\llLm > A^ 7o (l + 7o)7^ 7 o>Wp)<Xp) (8) 

where Yi m {p) is the spherical harmonic function and < JM\llLm > is the Clebsch-Gordan coeffi- 
cient. Substituting Eq.(|7|) and (|8]) in Eq.(|j), one derives the equations for the scalar wavefunction 

Hp) 

(i)S=0 

(M — Ei(p) — E 2 (p))gi(p)<j)(p) 
Ei(p)E 2 (p) + mim 2 + p 2 



AE 1 (p)E 2 (p) 

Ei(p)m 2 + E 2 (p)mi 
m (p) E 2 (p) 



d 3 k(G s (p - k) - 4Gy (p - k))gi(k)P L (cose)<f>{k) 



d A k(G s (p- k) + 2G v {p - k))g 2 (k)P L (cosO)0(k) 



+ E ^(p]e 2 (p) J d3kG s(p-k)p-kg 3 (k)P L (cosQ)m 

+ ^{p~)E 2 {p) I d ^ Gs{ -P~ $ ~ 2G v(P- *))P- kg 4 (k)P L (cose)m (9) 

where 

(Si(p) +m 1 )(E 2 (p) + m 2 )+p 2 



9i(j>) 
92(p) 



4£ 1 (p)£ 2 (p) 

{Ei(p) +m 1 )(E 2 (p) + m 2 ) - p 2 
AE x {p)E 2 {p) 
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(ii)S=l 



where 



Ei (p) + m 1 + £ 2 (ff) + m 2 
93[P) m( P )E 2 (p) 

_ Ei(p) + mi - E 2 (p) - m 2 
9dP) ~ ±Ei{p)E 2 (p) 

Ei{p) = yjp* + ml 
E 2 {p) = ^jp + ml 



{M-Ei{p)-E 2 {p))h{p)^p) 
= 4Ei(p)E 2 (p) { I d3k tt 2G v(P~ *) " °sip- k))fi(k)(mi + ma) 

G s (p - k)f 2 (k)(Ei{p) + £ 2 (p)]P L ( COS e)<£(A;) 
+ [ | d 3 fc(4Gy(p - fc) + k))h{k){Ei{p)E 2 {p) - mim 2 +f) 

+ (G s (p - fc) - 2Gy(p - k))f 7 {k)(miE 2 - m 2 Ei)]Pj(cose)^(k) 

+ | d 3 k[(2G v (p- k) - G s {p- k))h(k)(mi + m 2 ) 

- G s (p- k)h(k)(Ei + E 2 ))p- kPj(cosQ)^(k)}, (10) 



hip) = 4El{p ) E2[p) Wirt + Mi)(E 2 (p) + M 2 ) +/) 

/ 2 (P) = i^^ ((^(P) + ™0(3»(p) + - j?) 
2(£ 1 (p)+m 1 ) 



Hp) = U(p) = 
hip) = -kip) 



4Ei(p)E 2 (p) 
2 

~4Ei{p)E 2 (p) 



flip) = 4^)W ^ l(P)+mi ~^ (P) ~ m2) 

hip) = ^^¥) {Elip) + mi + E2{p) + m2) - 



The normalization condition / d?pTr{$> + (p)$>(p)} = -^p- for the BS wavefunction (f>(p) leads to 

To the leading order in the nonrelativistic limit, Eqs.(|7]) and @ are just the ordinary nonrelativistic 
Schrodinger equation for orbital angular momentum L with simply a spin-independent linear plus 
Coulomb potential. Solving equation (^) or (|T0|), we can get the spectra and wavefunctions for 
any given angular-momentum state 2S+1 Lj of heavy mesons. With these wavefunctions we can 
calculate hadronic matrix elements of the processes involving corresponding states, and the rela- 
tivistic corrections due to interquark dynamics are included automatically in them. This approach 
is different from convential ones which start from Schrodinger equation with all relativistic effects 
considered perturbatively. 



3 Decay rates of T(x c o — > 77) and r (Xc2 — > 77) 

Electromagnetic decays of XcO an d Xc2 proceed via the annihilation of cc to two photons. 
Here only electromagnetic interactions are considered, and color-octet components which contribute 
dominantly in hadronic decays of P-wave quarkonium do not contribute to electromagnetic decay 
widths, because final states are the photons which can not be produced via the annihilation of 
color-octet QQ pair. So two photonic decays of Xcj for J = 0, 2 can be well expressed in the BS 
formalism and relativistic corrections are incorporated systematicly in the decay rates. In the BS 
formalism the annihilation matrix elements can be written as follows 

< 0\QIQ\P >= J d 4 pTr[I(p,P) X p(p)}, (12) 

where I(p,P) is the interaction vertex of the QQ with other fields (e.g., the photons or gluons) 
which, in general, may also depend on the variable q° (the time-component of the relative momen- 
tum). If I(p,P) is independent of q° (e.g., if quarks are on their mass-shells in the annihilation), 
the equation can be written as 

< 0\QIQ\P >= J d 3 pTr[I(p,P)<S> P (p)}, (13) 
For process XcO ~^ 77 or Xc2 77 with the momenta and polarizations of photons k±,ei 
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and &2 , (-2 , the decay amplitude can be written as 

T =< 0\cT^c\ X cj > e^2 (14) 
for J=0,2, where ^1(^2) is the charm quark (antiquark) momentum, and 

i> = (?\iv-, — \ 7m + Tm ii — t iv] 

pi — fil — m fi — p2 — m 



Since P\+ P2 = M, as usual we take 



n n M , , 

Pl=P° = y ( 15 ) 



Therefore, the amplitude T becomes independent of p°. In terms of T, the decay rates can be 
written as 

nxcj -77) = ^EE / m a <« ( 16 ) 

spin polar 

for J = 0, 2, where the factor 1/2! is needed because N\ same graphs appear for N-photon final 
states. The photon polariztion is summed over in the Feynman gauge, 

X; ^(fc 1 )e w (fci) = -/r, 

Substituting BS wavefunction (Jsj) into (16), after summing over final states and averaging over 
initial states, we get 

r(XcO -» 77) = 24e^a 2 ( Cl + 3c 2 + 2c 3 ) 2 (17) 

12en« 2 n n 
r(x c2 ^77) = ^^(c 2 -2c lC3 + 7c|) (18) 



where 



J (p- k) 2 + m 2 2 2 

3^(p • k) _ 5(p- k) 4 4>(p) 
[ 4 2 4 Jl p 

C2 = [d*p—^ 2 — - 2 {^\?-(p-m+\\f-(P-m 2 }^ 

J [p — k) A + m z / 4 p 

c 3 = — -{ ~ E2 - mE if - (p- m + -\f - (p. km 



(p-k) 2 + m 2 2 4 p 
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In the nonrelativistic limit, (|17|) and (|i~8|) reduce to 



r(xco - 77) = — S-l / rfH»l 2 



32e 4 a 2 /■ 
r(x c2 - 77) = -g^T-l / ^ XC2 (P)| 2 



5m' 

Using the Fourier transformation of wavefunctios 

/ d 3 pp<t> XcJ (p) = ^=KcA°) 

we derive the well known results in coordinate space, which is consistent with that given in Q 

27e 4 a 2 

r(x c0 -77) = ^jH*U°)l 2 ( 19 ) 
r(xc-77) = ^^l^ c2 (°)l 2 ' ( 2 °) 

where R' j(0) is the derivative of radial wavefunction at the origin, and in the nonrelativistic limit, 
-^Xco(O) = ^Xc2^)' due ^° ^ e heavy quark spin symmetry. 

Recently, in the framework of NRQCD the factorization formulas for the long distance and 
short distance effects were found to involve a double expansion in the quark relative velocity v 
and in the QCD coupling constant a s [13, 14]. To next-to- leading order in both v 2 and a s , as an 



approximation, we may write 

r(XcO - 77) = 244a 2 ( Cl + 3c 2 + 2c 3 ) 2 [l + ^(y - f )] (21) 

r( X c2 - 77) = l -^f-(c 2 - 2c lC3 + 7c 2 )(l - ^-j ) (22) 

where we have used QCD radiative corrections given in [|l5|]. We must emphasize that above 
factorization formula are correct only to next-to-leading order in v 2 and a s . If higher order effects 
are involved, the decay widths can not be factored into a integral of wavefunction and a coefficient 
that can be written as a series of a s . NRQCD has applied a more general factorization formula for 
quarkonium decay rates, which will be discussed in detail later. 
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For the heavy quarkonium cc systems, m\ = m<2 = m c , Eqs. (B) and (|10|) become much 
simpler. We take the following parameters which appear in potential (|5|), 

m c = 1.5Gev, X = 0.23Gev 2 , A QC d = 0.18Gev, 

a = OMGev, a = e = 2.7183. 

With these values the mass spectrum of charmonium are found to fit the data well. In Fig.l and 
Fig. 2 the solved scalar wavefunctions both in momentum and coordinate space for P-wave triplet 
Xcj states are shown and we can see explicitly the differences between wave functions for J = 0, 2 
but they are same in the nonrelativistic limit. Substituting <^ Xc0 (p) and 4> Xc 2(p) m t° ( HI) an d (H)i 
we get 

T(xo -» 77) = 5.32keV, 
T(X2 -» 77) = 0A4keV 

their ratio is 

R = r j X0 ^ 77 j = i2.l. (23) 
T(x2 -* 77) 



Our results are satisfactory, as compared with the Particle Data Group experimental values [16| 
r(XcO - > 77) = 5.6 ± 3.2keV, and T(xc2 77) = 0.32 ± O.lkeV. Here in above calculations the 
value of a s (m c ) in the QCD radiative correction factor in (21) and ( p2[ ) is chosen to be 0.29 Q, 



which is also consistent with our determination from the ratio of B(J/ip — > 3g) to B(J/ip — > e + e ) 
1 • 

Moreover, in order to see the sensitivity of the decay widths to the parameters, especially 
the charm quark mass, we use other two sets of parameters 

m c = lAGev, A = 0.24Get> 2 , 

m c = 1.6Gev, A = 0.22Get> 2 , 

with other parameters keeping unchanged (the heavy quarkonia mass spectra are not sensitive to 
a, a for a < 0.06GeV), By the same procedure, we obtain 

T(xo -» 77) = 5.82(4.85)/cey, 

T(X2 -» 77) = 0.50(0.39)A;ey, 
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and the ratio 



R = ~* 77 j = 11.8(12.5) (24) 



for m c = 1.4(1.6)GeV. 



We find the widths are decreased with the decreasing of A This is obvious since the wave- 
function in coordinate space will become broader when the slope of linear potential is decreased 
and the corresponding wavefunction in momentum space will become narrower, so the effective 
couplings decay become smaller. It is interesting to note that the ratio of two photonic decay width 
of XcO and Xc2 is almost unchanged and is insensitive to the choice of parameters. 

Finally we discuss the relation between our approach and the NRQCD theory. Recently, a 
general factorization formula which is based on nonrelativistic QCD (NRQCD) has been developed 
for studying the inclusive cross sections of production and decay of heavy quarkonium. In this 
formalism the quarkonium decay rates can be written as a sum of a set of matrix elements to any 
given order in v 2 , with each matrix element multiplied by a coefficient which can be calculated in 
perturbative QCD. This approach has been proved successful in the application of some processes 



involving heavy quarkonium [17, [Eg]. In NRQCD, the electromagnetic decay rates of XcO and X02 
to next-to-leading order in v 2 can be written as 



tv n 2Imf EM ( 3 Po) f3 „s, 2Img EM ( 3 Po) ^ , 3 „», 
TUo -> 77) = 4 < Xco\Oem{ PoHXco > + s < Xco\Oem{ PoHXco > 



(25) 



tv ^ 2/m/ EM ( 3 p 2 ) ^ l/n ,3 dm ^ 2Img E M{ 3 P2) . un , 3e >m 

T(X2 -> 77) = 4 < Xc2 \Oem{ P2) \Xc2 > + s < Xc2\Oem{ P2)\Xc2 > 



where 



(26) 



Oem( 3 Po) = ^ + (-^D)-ax|0><0|x + (-^D)-^ 

Oem( 3 P2) = B)(V i) x|0 >< 0|x+(-~ 

Gem( 3 P ) = ^[^ + (-^D) 2 (-^D)-<?x|0><0|x + (-^D)-^ + / i .c] 

Gem( 3 Po) = \[^ + {-\ D) 2 (-^ Dfa^xlO >< 0|x + (-| B) W + h.c] (27) 
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where D is the space component of covariant derivate D^, tp and \ are two component operators 
of quark and antiquark respectively. If identifying the quark operator expectation values with the 
deriatives of wavefunctions at the origin, the decay widths can be written as 



r(xo 77) 



+ 1Mm Zf n) Re(R^ m m (28) 



■nrrr 



9Imf E M( 3 P2) lT} , , m|2 . 15Img E M( 3 P2) D , D ( 3 )/ m w /„« / onA 
T(X2^77) = —4 I^WI + =e (0)^(0)) (29) 



7rm u 



In comparision, we take the on-shell condition p\ = p\ = E instead of (|T^) and expand 
the annihilation amplitudes in (|l7|) and ( |l8|) in terms of p 2 /m 2 . The leading order contribution 
of P-wave decay comes from terms linear in p and in order to take into account their higher order 
effects we retain those terms proportional to the third powers of p then have 



3e%a 2 r „ # 
r(xo-77 =-Vl d 3 pp(l 1 



rrr 



, Qs ,7T 28., 

f , m2 ^)l 2 [i + -( T - ¥ )] 



(30) 



■ie%a 2 r „ „ a« 16 

T(X2 - 77) = -g^H y rfVfeb)| 2 (l - ^y)) 



(31) 



where the stanard Schrodinger wavefunction (with relativistic corrections) 4>Sch(p) is related to 4>(P) 
through the normalization condition (11) 



4>sch{p) = -^=( m ~l E )(t>(p), 
VM E 



(2vr) 3 y dpp^sMl 2 = 1. 



Using formula 



d 3 pp4>Sch(p) = 3R' Sch (0), 



d 3 pp 3 <l>Sch(p) = 541(°)> 



the expressions of fl30| ) and fl3l|) are transfered into coordinate space and comparing with that 



derived from NRQCD (|28| ) and (2£), we can easily determine the coefficients 



Imf EM ( 3 P Q ) 



o 4 2n , a s (K 2 28 

37re?)Q! 1 H ( ) 

vr 3 9 ;J 
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„ 4ire%a 2 a< , 16 

Imf EM eP 2 ) = — 

Img EM ( 3 P 2 ) = (32) 
Here we only consider the QCD radiative corrections to leading order coefficients ImfEM^Po) an d 



^™/bm( 3 ^) which are equal to the results derived in [13|. More, we have determined two new 
coefficents, i.e, the second one and fourth one in (^). These two matrix elements have a suppression 
factor of v 2 so we need not to take into account higher order corrections to their coefficients any 



more. 



4 Summary and Discussion 



In this paper we provide an estimate for the photonic decays of P-wave charmonium. It is clear 
from above calculations that comparing with nonrelativistic results the relativistic effects enhance 
the ratio R substantially. We know that there are two sources of relativistic corrections: l)the 
correction due to relativistic kinematics which appears explicitly in the decay amplitudes; 2)the 
correction due to inter-quark dynamical effects (e.g. the well known Breit-Fermi interactions), which 
mainly causes the correction to the bound state wave functions. From the expressions (^) and (|3~I|) 

of decay rates which have been expanded to the firdt order of one might expect that the ratio 
R would become smaller after taking relativistic corrections into account, because the coefficient of 
the term ^ in (|30| ) is smaller than that in ([n]). However, mainly due to the attractive spin-orbital 

force induced by one gluon exchange for the ++ meson, the xo wavefunction becomes narrower 
than \2 wavefunction in coordinate space, and therefore the derivative of wavefunction at the orgin 
become larger for xo than that for X2- As a result, the dynamic relativistic effect on R is in the 
opposite direction and can be even larger. The overall relativistic correction to R is found to be 
positive, and our result is in agreement with the E760 data and disagree with the values mesured 
by CLEO and TPC2. 

Our expressions for the decay widths are identical with that derived from the rigorous 
factorization formula to next-to- leading order in v 2 and in a s . Moreover we have determined 
two new coefficents in the nonperturbative matrix elements for these decay widths. For a more 
accurate estimate, higher order corrections both in v 2 and in a s should be taken into account. For 
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electromagnetic decays, in general we can estimate them within the \QQ > sector and avoid the 
difficult problem due to the effects of high Fock states such as \QQg > and \QQgg >■ But we must 
notice that if higher order matrix elements are included, the decay widths can not be factored in the 
way like (^) and (p^) because the higher order coefficients are different for each nonperturbative 
factor. 

We have solved the BS equation for the bound-state wave functions with QCD inspired 
interquark potentials (linear confinement potential plus one gluon exchange potential) as the BS 
kernel. With some popular parameters for the potentials we obtained the wave functions and used 
them to calculate the decay widths. From (^) and ( p4| ) it can be seen that three different wave 
functions lead to somewhat different photonic decay widths but give very close values for R. This 
might indicate that our estimate of R is insensitive to the quark mass and potential parameters, and 
therefore could be a rather reliable result, despite the uncertainty in the estimate of the dynamical 
relativistic effects. We hope the lattice simulations will give more reliable estimates for these decays 
within the framework of NRQCD, and can be compared with our results. 
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(dashed line) and Xc2 (dashed and dot line) by solving BS equations with m c = l.5GeV 

Fig.2 Wave functions Rs c h(%) (normalized in coordinate space) of XcO (solid line), Xd 
(dashed line) and Xc2 (dashed and dot line) by solving BS equations with m c = 1.5GeV 
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